Projecting Fermion Pair Condensates into Molecular Condensates 
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We offer strong evidence that the recent observations by M. Greiner, C. Regal, and D. JinpjJ and 
by MIT group are signatures of a fermion superfiuid in the strongly interacting regime made up 
of large fermion pairs. Our conclusions are based on calculations using crossover theory for different 
potentials including those with the characteristics of two-channel models. Our results demonstrate 
clearly universality near resonance. The T c predicted by crossover theory is a perfect match with 
the observed boundary of vanishing condensate fraction with no adjustable parameters. 



In an exciting recent paper, M. Greiner, C. Regal, 
and D. Jin reported evidence for condensation of fermion 
pairs of 40 K near Feshbach resonance 0. Similar results 
in 6 Li were reported later by Ketterle's group at MIT 
but with a different interpretation^- Feshbach reso- 
nance is produced by Zeeman shifting the energy of a 
bound state in closed channel to zero energy. When this 
energy (linear in magnetic field B and usually referred 
to as "detuning" 5) vanishes, a low energy fermion pair 
in the open channel strongly resonates with the bound 
state through hyperfine interaction, leading to consider- 
able scattering. The regions S > and S < (which has 
negative and positive scattering length a s respectively) 
are referred to the BCS and BEC side of the resonance 
because the ground state deep inside these regions are 
BCS superfiuid and Bose-Einstein condensate (BEC) of 
molecules respectively. The intriguing region is near res- 
onance (S — > 0) where the system becomes strongly in- 
teracting. The experiments in ref.pj and Q explore this 
region. 

In ref. lj , an equilibrium state was prepared at an ini- 
tial field Bi, which was quickly changed to Bf deeper in 
the BEC side. It is found that when Bi is sufficiently close 
to resonance, a molecular condensate will emerge at Bf 
such that the fraction of condensed molecules (N /N m ) is 
independent of the (fast) sweep rateQ. The JILA group 
considers this as evidence for condensation of fermion 
pairs in the initial state. Their picture is the BCS-BEC 
crossover theory, which describes the BCS to BEC evo- 
lution as a smooth reduction in size of the condensed 
pair 0,0,0. The immediate question is that near reso- 
nance the size of the pair should be of interparticle spac- 
ing, then the direct overlap between these pairs and the 
tightly bound molecules at Bf should be very small. In 
that case, the final molecular condensate must be gener- 
ated by the fast sweeping field, and the relation between 
the momentum distributions in the final state and the 
initial state is not immediately clear. 

In contrast, the MIT group thinks that their similar 
results are due to a BEC of tightly bound molecules in 
closed channel persisting onto the BCS side. They argue 
that even though these molecules are unstable in vac- 
uum on that side, they can be stabilized by the Fermi 
sea. Somewhat earlier, Falco and Stoof (FS)@ had 



also suggested that the molecular fraction N /N m ob- 
served at Bf in ref. 1] is given by the molecular com- 
ponent Z{8) in the condensed fermion pair in the ini- 
tial state, | pair) = \JZ{5)\ closed) + \J\ — Z(<5)|open), 
where | closed) and |open) are closed channel and open 
channel fermion pairs. Using a phcnomcnological reso- 
nance model, FS calculated the boundary of vanishing Z 
and found reasonable match with the boundary of van- 
ishing N /N m observed in ref.0. FS argue that even 
though when Z vanishes, the system still has pairing in 
open channel, and the observed boundary of vanishing 
No/N m is a sharp crossover from molecular rich pair to 
fermion pair, (see fig-1 in ref.Q). This "persisting molec- 
ular BEC" pictures^, El are very different from the usual 
crossover picture, which does not predict a sharp division 
between molecule rich and pure fermion pairs. 

Adding to this confusion is the question of whether the 
many-body physics near a Feshbach resonance can only 
be described by "two channel" (or "resonance") mod- 
els which incorporate open and closed channel physics 
specifically 0, or can be described by simpler "single 
channel" models which consist of two types of fermions 
interacting with a tunable potential|3|, |J|. Should these 
models have very different properties at resonance, and 
if the physics near resonance is universal as generally be- 
lieved, it would imply that there are two distinct univer- 
sality classes, and that only the two channel universality 
class gives the right physics of Feshbach resonance. The 
other possibility is that two channel models are not uni- 
versal, and that the notion of universality inferred from 
current experiments[2ll3| is misconstrued. 

In this paper, we point out that these confusions can 
be eliminated by deriving the precise expression for the 
observed molecular fraction, calculating the properties of 
the system explicitly using crossover theory, and examin- 
ing the nature of the broken symmetry. We first perform 
T = calculations using crossover theory [2l |4| for three 
different types of potential: Vi : a square well, Vii'. a well 
with a high barrier, and Vm a delta-function potential. 
Since Vu can accommodate a very long lived quasi-bound 
state, it has the major characteristics of a two channel 
(or "resonance" ) model. We shall show that all these po- 
tentials give rise to identical properties (such as chemical 
potential, energy gap, and coherence factors) near res- 
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onance, demonstrating universality in this regime. The 
fact that Vii carries all the characteristics of the two- 
channel model suggests both types of models may have 
identical properties near resonance. To our surprise, our 
calculations show that there is non-trivial overlap be- 
tween the pairs at resonance (denoted as i = 0, where 
x = (kpa s ) arid kp is the Fermi wavevector) and those 
at x as large as 12, which corresponds to the the final field 
Bf in ref.jlj. This eliminates a serious concern formerly 
expressed by one of us (TLH) 9] about the validity of pro- 
jection argument in ref.Q]. To the extent that direct pro- 
jection describes the experimental process, one can show 
that No /N m is directly proportional to the superfluid or- 
der parameter, and the boundary for vanishing No/N m 
in T — B plane observed in ref . [jj, |2| is precisely the T c 
between superfluid and normal gas, not a crossover be- 
tween different types of pairs. Since Vj to Vm have iden- 
tical properties near resonance, one can use the simplest 
potential (Vm) to study the finite temperature proper- 
ties in this region. Such study had been performed in 
ref.Q by studying Gaussian fluctuations about the mean 
field state. Remarkably, the phase boundary predicted in 
ref.0 matches very well with the observed boundary for 
vanishing condensate fraction (No/N m = 0)0,0 with no 
adjustable parameters. 

(A) The fraction of condensed molecules: Con- 
sider a fermion pair with total momentum q, Dt(x) — 

Ek,^/k,^( 3; )4 + q/2, a tt -k + q/2,/3/ 2 ' where °k,a creates 
a fermion with momentum k and hyperfine spin a, 
x = (kpa s ) describes the distance from resonance, and 
/k ad{ x ) i s t ne Fourier transform of the pair wavefunction 
f a0 (r;x) = Q~ 1/2 Y,k elkr fk, a p( x ) which depends on x; 
O is the volume of the system. A condensate of N zero 
momentum pairs is |x) = Af D]£[ (x)\vac) . It reduces to 
a weak coupling BCS superfluid and a molecular conden- 
sate as x « — 1 and x » 1 respectively^^!]. In grand 
canonical description, this pairing state becomes 

l*0»0> =M II f>(^) + Wk,a/s(a5)4,a a -k,/9) l V&C > (*) 
k,a/3 

where Af is the normalization constant, u and v are 
related to the Fourier transform of f a /3(f) as fk,a/3 = 
C v Kap/ u k: an d C is a normalization constant. 

Let x and x be the distances of the initial and fi- 
nal state from resonance, and x be on the BEC side. 
When the system is jumped from x a to x, the number of 
condensed and un-condensed molecules emerge at x are 
N = (Dl(x)D (x)) Xo and N ex = E^oW^Aj 
respectively. We then have from eq.jlj, 



N = \(D (x))f 
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N K 



; Ek, q l/k,a/3(a;)| 2 n- q , Q n q _2k,/3/2, where ^ a(3 (x ) = 
( a k.Q a -k,/3)^o: "k = ( a k .a a k,a)zo> and we nave ignored 



in eq.© a term down by a factor of iV _1 . Eq.f$) shows 
that Nq is the overlap of the initial order parameter \t a p 
with the final pair wavefunction. 

Let us first discuss the single channel case because it 
is illuminating. Denoting the two spin states as j and J,, 
the s-wave pairing function is f a p(r) = f(j)io~ v a p, and 



No 

N m 



IEk/kW^ k | 2 

I Ek /k (x)n\ 2 + Ek, q l/k(x)| 2 n q n q _ 2k 
. = N + N ex , *(k)* = (4t°-kl)».. 



(3) 



where N„ 

( a ic | a k = ( a k i a k i)x - Since 'J'k vanishes at T = T c , 
the curve Nq /N m — in the T — B plane is therefore the 
boundary of superfluid to normal transition. 

(B) The Crossover Picture: The basic assumption 
of crossover theory is that the ground state near reso- 
nance is still given by eq.JQ). We have studied eq.JSJ us- 
ing the single channel crossover theory for: (I) a square 
well, Vi(r) = —\U \ (or 0) for r < r D (or r > r Q )\ 
(II) a well plus a barrier, Vn(r) = —\U„\ for r < r a . 
Vn(r) = \Ui\ for r Q < r < n, and Vzr(r) = for 
r > ri; (III) a contact potential Vm(r) — gS(r), where 
g = 2-kT% 2 a s / 'M , where M is the mass of the fermion. The 
Hamiltonian is H = E k ,a e k4,a a k,c* + Ek,k', q V ( k ~ 

k/ ) a q /2+k,T a q/2-k,i a q/2-kM a q/2+k',i' e k = ^ k 2 /2M. 
The range of the potential of Vi and Vn (r and n) 
are taken to be of atomic size, much smaller than inter- 
particle spacing, i.e. k F r a ,k F ri << 1. By varying \U a \, a 
bound state can be peeled off from the continuum, caus- 
ing a s to jump from — oo to +oo. The main difference 
between Vj and V/j is that the latter can accommodate 
a quasi-bound state of atomic size, similar to the feature 
of the two channel or the resonance model. 

For all cases, itk and v are given by the well known 
expressions |u k | 2 = (E k + £ k )/(2£ k ), K| 2 = (E k - 

M; #k 



4)/(2^ k ), where £ k = e v 



and 



the energy gap A k at T = is given by the gap equation 
A k = - Ek' V ( k - k ') A w/( 2E k')- The corresponding 
expressions for \I> k and nk are ^ k = u k v k — A k /(2£'k), 
rik = f k . The chemical potential /i is determined by the 
number density n as n = f2 _1 Ek n k~ = 0' I 1 )- 

(C) Results: We have solved the gap equation nu- 
merically and have inverted the relation n = n([x) so that 
quantities like A k , u k , v k , ^ k , n k are functions of the pa- 
rameters of the potential (denoted collectively as {!*}) 
and density n. Our results are displayed in figure 1 to 
9. In these figures, we use the dimensionaless parameter 
x = (kpOg)^ 1 instead of potential parameters {Yi}. The 
results for potential Vj ,Vn, and Vm will be colored in 
black, red, and green respectively. For both Vj and Vn, 
we take kpr D = 0.05. For Vu, we choose r\/r = 1.2 and 
U x = 5|C7 (min) |, where |C/^ mm) | = h 2 ir 2 / (4Mr 2 ) is the 
value for a square well to produce a bound state. 

(CI) At, A , and /i: In figure la, we have plotted 
Ak/A versus kr Q for x = 2, 0, —2, denoted as dashed, 
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FIG. 1: Fig.la to lc: Black, red, and green represent Vi,Vii,Vni, x — (kFds) -1 , and kpr — 0.05. Dashed, solid, and dotted 
lines denote x — 2,0, —2. Within the range —2 < x < 2, there are little changes in A(fc)/A(0), although A(0) and fj, undergo 
significant changes (see fig. lb), leading to very different ^(k) in fig. lc. 



solid, and dotted lines respectively. The inset shows that 
Ak/Ao for different potentials have similar values at low 
momenta (even up to a few fc^'s) but differ significantly 
at higher momenta, though it vanishes for realistic po- 
tentials Vi and Viz- Within the range — 2 < x < 2, this 
ratio is essentially unchanged. The value of A(0) and fi, 
however, changes a lot as x varies from -2 to 2, as shown 
in fig. lb. These variations lead to the large differences 
in \I/ k as one move across the resonance. It is clear from 
fig. lb and lc that universal behavior sets in around res- 
onance (x — 0) since different potentials give identical 
results. For all potentials, we find that nt is accurately 
described by that of Vnj. Its figure is given in ref.Q. 

(C2) Pair wavefunction /: To display /, we plot 
■J 4ir/kFrf(r) versus kpr in fig. 2a to 2c. One sees that 
all potentials have essentially the same pair wavefunc- 
tions unless one goes deeper into the BEC side. The 
function / on the BEC side for all potentials have essen- 
tially the Feshbach form oc e~ T l as jr. Until a s reduces 
to the atomic size r a , f has a non-trivial fraction be- 
yond r a . It can therefore have significant overlap with 
the pair wavefunction near resonance. In fact, the over- 
lap between / at x = and that at x — 12 (which is the 
final state in ref. 1]) is as big as 18%. Thus, the direct 
projection picture suggested in ref. Q] (formerly criticized 
by one of us 9]) is in fact reasonable. In the case of Vn, 
we note that while a molecular component of range r Q 
is visible on the BCS side near resonance, it essentially 
vanishes for x < —2. However, the T = condensate 
fraction No/N m is still very large for initial state x = —3 
which has no molecular component. This shows that the 
molecular component is unimportant in contributing to 
the projected condensate fraction No/N m . 

(C3) The condensate fraction No/N m , and the 
superfluid-normal phase boundary: In figure 3a, we 
show the T = molecular fraction No/N m at x = 12 
as a function of the initial state x . The data of ref. [J 
and are represented by circles and square respectively. 
The former is obtained at T/Tp = 0.05. The temper- 
ature of the latter is uncertain. Since the decrease of 
No/N m on the BEC side appears to be related to three 



body effects [l(, we apply equilibrium theories only to the 
BCS side {x < 0). The large difference in scale between 
the results at T = and at T/Tp = 0.05 is anticipated 
since TV is strongly suppressed in the strongly interacting 
regime [3,0, which will reflect in a strong suppression of 
No/N m . The full effect of the crossover theory is best 
illustrated by showing T c in the T — x plane. In fig. 3b 
and 3c, we have superposed the T c predicted in ref.Q 
on the contour plots for No/N m in the JILA^(| and the 
MIT 2] experiments. (Because of the universality near 
resonance demonstrated in fig.l and 2, we can use the 
finite temperature results for Vjjj |J| to compare with ex- 
periments.) The matches are remarkable, in view that 
there are no adjustable parameters. 

In fig. 3a, we also show that as the final state moves 
deeper into the BEC side, the boundary of vanishing 
No/N m moves to the opposite, i.e. BCS side. This is 
because while both No and N ex decrease as the final 
state becomes more tightly bound, the decrease of N ex is 
faster than that of Nq so that No/N m actually increases. 
This fact can be demonstrated numerically, but can be 
established analytically if we take f(k) — (k 2 + A 2 ) -1 
and study the change of No/N ex as a function of A 2 . 
As A — > oo j No/N m reaches the asymptotic curve T: 
No/N m = A/(A + B), where A = IE k *kl 2 > and 
B = Eqk n k n k-qi a prediction that can be verified ex- 
perimentally. Since jumping to fields where molecular 
size are comparable r a (such as a deeply bound state) 
means sampling short distance behaviors of the order pa- 
rameter (see fig.la), T will not be universal, as illustrated 
in the dashed and dotted curve in fig. 3a. 

(D) Two channel models: For Feshbach reso- 
nances, in Section (A) is f a p{r)=\fZf^ (r)x^(r) 

+Vl — Z f^°\r)x^°p{r), where / and x are normal- 
ized orbital and spin functions, "c" and "o" denote 
close and open channel, and and x are or- 

thogonal. The zero momentum pair is -Dq =0 = 

VZbl =0 + VT^E k / k 0) xa«L«l k ,^ where b^ = 
J2k / k C ^XQ/3 a kQ a L k p is the closed channel boson in 
the resonance model|]|. From eq.J2J, we have 
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FIG. 2: Fig. 2a to 2c: Pair wavefunctions for different values of x, a — ^/47r/fcF. For Vu, f has a peak at r < r a , yet most of 
the weight for all potentials is outside the atomic range r . At resonance the size of / is given by fc^ 1 . 




FIG. 3: Fig^a: The T = value of No/N m at x — 12 as a function of initial state x is represented by the solid line. The 
data of ref.pj and are denoted as circles and squares. As the final state x moves toward the BEC side, the curves moves 
toward the BCS side, reaching a fixed asymptote, which depends on the microscopic nature of A(fc) at large momenta (see 
fig. la). The asymptotes of Vi and Vu are represented by the dotted and dashed curve resp. In fig. 3b and 3c, we superpose 
the phase boundary calculated in ref.(3J (red triangle) on the data of ref.[Tfl] and 0. In fig. 3b, light color means large No/N m - 
The numeral correspondence of the color code is given in ref.|lOl|. Dark blue corresponds to No/N m = 0. In fig. 3c, we have 
used the approximate formula x = (832 — B)/(B — 640) to convert magnetic field into x — ffcFa«')~ 1 [Tll|. 



*o = iV^Ek Jf*r + VT^Ek/i 0) *ri a . w ^e 
^,(c,o) _ Trx^ c '°^^ k /2 are the components of the order 
parameter ^k,a0 in closed and open channel. If the fi- 
nal state is a tightly bound molecule (Z = 1), then the 
vanishing of N a is given by the vanishing of ^^ c \ How- 
ever, to infer from a vanishing of a possible non-zero 
\?(% (like the vanishing Z in / means the open channel 
component (1 — Z) — ► 1), is incorrect. This is because 

(c) 

^k,a/3 has no normalization constraint, and both "f^ 

and ^ are dynamically connected through the coupling 

term frq = o a ka a -k / 3X/3° c ^ + ^- c - i n the Hamiltonian. Any 
non-vanishing pairing order in one channel will generate 
a non-vanishing order in the other. 

We have also learned from D. Jin that for 40 K there 
is a shallow bound state (C) in the open channel. Thus, 
as the close channel bound state (23) is Zeeman shifted 
down, it will cross the zero energy state (A) and then 
the bound state (C) in the open channel. One can then 
go from A to 23 to C using a magnetic field sweep, so 
that a pair of fermions in the open channel can be swept 
into a bound state in the same channel. In that case, the 
system can be described by single channel models. 
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